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Biorthogonal  Refinable  Spline  Functions 


Tim  N.  T.  Goodman 


Abstract.  We  give  a construction  for  refinable  spline  functions  of  degree 
n with  compact  support  and  simple  knots  in  \7L  which  are  biorthogonal 
to  uniform  B-splines  of  degree  n with  simple  knots  at  \TL. 


§1.  Introduction 

A function  is  refinable  if  it  is  a linear  combination  of  dilates  of  integer  translates 
of  itself.  Such  functions  are  central  to  multiresolution  methods,  in  particular 
in  the  construction  of  wavelets.  In  general,  refinable  functions  can  be  defined 
implicitly  from  the  refinement  equation  which  they  satisfy,  but  explicit  con- 
structions of  refinable  functions  are  restricted  mainly  to  spline  functions,  i.e. 
piecewise  polynomials.  If  we  require  the  natural  condition  that  the  integer 
translates  of  the  univariate  refinable  spline  function  <j>  with  compact  support 
form  a Riesz  basis,  then  cf>  can  only  be  a uniform  B-spline  with  simple  knots 
[4].  However  there  is  more  flexibility  if  we  replace  the  single  function  <f>  by 
a refinable  vector  of  spline  functions  . . . , <j>r).  For  a survey  on  refinable 
spline  functions,  see  [2], 

In  multiresolution  methods,  orthogonality  plays  an  important  role.  In  [1], 
constructions  are  given  for  refinable  functions  whose  integer  translates  are 
biorthogonal  to  a given  refinable  function  <j),  in  particular  when  <j>  is  a uniform 
B-spline  with  simple  knots.  However  these  dual  functions  are  not  defined 
explicitly.  We  give,  in  Section  3,  constructions  for  refinable  spline  functions 
of  compact  support  which  are  biorthogonal  to  uniform  B-splines  with  simple 
knots.  This  requires  refinable  vectors  of  three  functions:  the  uniform  B-splines 
have  knots  in  while  the  dual  functions  have  the  same  degree  and  simple 
knots  in  |5Z.  The  construction  is  based  on  a general  result  in  Section  2 giving 
necessary  and  sufficient  conditions  for  biorthogonality  of  certain  vectors  of 
(not  necessarily  refinable)  functions  in  terms  of  a Grammian  matrix. 
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§2.  Biorthogonal  Basic  Sets 

Let  <p be  compactly  supported  real- valued  functions  in  L2(1R).  We 
say  . . . ,<f>r}  is  a basic  set  for  a space  V if  V comprises  all  real,  finite, 
linear  combinations  of  integer  translates  of  (pi, . . . , 4>r. 

We  say  basic  sets  {(pi, . . . ,(pr}  and  {ip i, . . . ,ipr}  are  biorthogonal  (or  the 
basic  set  {ipi, . . . ,ipr}  is  dual  to  {(pi,  . . . ,<pr})  if 

f°° 

I 4*i'Pj{‘  i)j  — 1,  • • • ,T,  k £ r7Zi. 

J — oo 

A basic  set  {(pi, . . . ,<pr}  is  said  to  be  stable  if  {</>*(.  — j)  : i = 1, . . . , r,j  £ 
7Zi}  forms  a Riesz  basis,  i.e.  for  some  A,  B > 0, 

r oo  -oo  r oo  r oo 

E <•«</  C E <■««•- a2  <bE  E 4 

i=l  j=— oo  00  1=1  j=—oo  i=l  j~  — oo 


for  any  a*j  € 1R,  2 = 1,...,  r,  j € ZZ. 

It  is  shown  in  [3]  that  {(pi, . . . ,<pr}  is  stable  if  and  only  if  for  each  u in 
1R,  there  are  integers  ki,...,kT  with 

det  ^4>i{u  + 2-7rfcj)j  ^ 0.  (1) 

We  shall  say  a matrix  M(z ) = [M(z)jj]|4=1  of  Laurent  polynomials  is 
invertible  if  it  has  an  inverse  which  is  a matrix  of  Laurent  polynomials,  i.e. 

det  M{z)  = azl,  some  a ^ 0,  l € ZZ. 

Lemma  1.  If  {<f> i, . . . , cpr}  is  a stable  basic  set  for  V,  then  { ipi , . . . , ips}  in  V 
also  forms  a stable  basic  set  for  V if  and  only  if  s = r and 

T OO 

'Pi  = £ £ - k)>  * = 1,  - - - » (2) 

j = l k~  — oo 


where  the  matrix  of  Laurent  polynomials 


A(z)  := 


£ Aim* 

.k=  — oo 


r 


i,j  = 1 


is  invertible. 

Before  proving  this  lemma,  it  will  be  useful  to  introduce  the  following 
vector  notation.  For  a basic  set  {(pi, . . . , </y},  we  let  <p  denote  the  column 
vector  {pi, . . . , cpr)T . Then  we  can  write  (2)  as 

OO 

k=—o o 


(3) 
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where  Ak  is  the  matrix  [Aij(k)}rl]=1 . Taking  Fourier  transforms  then  gives 


'>/’(«)  = A(z)<j>(u), 


(4) 


where  z = e_,“. 

Proof  of  Lemma  1:  Suppose  that  (2)  holds,  where  A(z)  is  invertible.  From 
(4)  we  have 

4>{u ) = A(z)~1ip(u). 

Since  A(z)_1  is  a matrix  of  Laurent  polynomials,  it  follows  that  <j) l,  ■ • • , <j>r 
are  finite  linear  combinations  of  integer  translates  of  ipi , . . . , ipr . Since  V 
comprises  all  finite  linear  combinations  of  integer  translates  of  <j>i, ...  ,<frr,  it 
follows  that  {-0i , . . . , ipr}  is  a basic  set  for  V. 

Also  for  any  u in  R,  we  may  choose  integers  ki, . . . , kT  so  that  (1)  holds, 
and  thus  from  (4), 

det[i/),(u  + 2'rrkj)}rij=1  = det  A(z)  det[^(u  + 27r/cj)]-J_1  ^ 0. 


So  {ipi, . . . , Vv}  is  stable. 

Conversely  suppose  that  {^i, . . . , ips}  is  a stable  basic  set  for  V.  Then 
there  exist  ansxr  matrix  A(z)  and  an  r x s matrix  B(z)  of  Laurent  polyno- 
mials such  that  for  z = e~,u, 


ip(u)  = A(z)(j>(u) 


and 

4>(u)  = B(z)ij}{u)  = B{z)A{z)4>{u). 

For  any  u € R we  may  choose  integers  k\,...,kr  so  that  (1)  holds.  Since 

\4>i{u  + 2nkj)]ij=1  = B(z)A(z)[4>i{u  + 2 Trkj)]Tij=1, 

it  follows  that  B(z)A(z)  = Ir,  the  r x r identity  matrix.  Similarly  V>(w)  = 
A(z)B(z)/ip(u)  and  since  {-0i, . . . , V's}  is  stable,  we  can  deduce  as  above  that 
A(z)B(z)  = Is.  Thus  s — r and  A(z)  is  invertible.  □ 

Theorem  2.  Suppose  that  {<f>  i, ...  ,(j>r}  and  {0i, . . . , ipr}  are  stable  basic  sets 
for  V and  W respectively.  For  k e 7Z,  i,j  = 1, . . . ,r,  we  define 

Mij(k)  :=  f0  - k),  (5) 

J — OO 

and  let  M denote  the  r x r matrix  of  Laurent  polynomials  given  by 

r 

(6) 

*Sj  = 1 


M(z) 


E 


Lk=- 
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Then  there  exist  biorthogonal  basic  sets  for  V and  W if  and  only  if  M is 
invertible.  Moreover  in  this  case,  for  any  stable  basic  set  for  V there  is  a 
unique  dual  stable  basic  set  for  W. 

Proof:  Let  {<pi, ...  ,<pr}  and  {ip i, . . . , ipr}  be  any  stable  basic  sets  for  V and 
W respectively.  Then  by  Lemma  2.1  we  have 

oo  oo 

0=  E *=  E (7) 

k~ — oo  k=  — oo 

where  Ak,  Bk  are  r x r matrices  such  that 

OO  OO 

A(z):=  £ Akzk,  B(z)  :=  £ Bkzk,  (8) 

fc=  — OO  k~  — OG 


are  invertible. 

For  k £ 7L,  i,  j = 1, . . . , r,  define 

Mn(k)  = r Mj{.  ~ k),  (9) 

J — OO 

and  let  M denote  the  r x r matrix  of  Laurent  polynomials  given  by 

r 

(10) 

».i=i 

Then  from  (2.5)-(2.10)  we  have  for  z = e~'u, 

M(z)  = A(z)M(z)B(z)* , (11) 

where  B(z)*  = B(zf  = B{z~})T . 

Now  by  (2.9)  and  (2.10),  4>  and  ip  are  biorthogonal  if  and  only  if  M{z ) = I. 
If  this  holds,  then  by  (11),  M must  be  invertible.  Conversely,  if  M is  invertible, 
then  for  any  choice  of  A{z)  we  can  define  B(z)  uniquely  by 

B{z)  = {M{z)-1A{z)-1y 

so  that  (11)  holds  with  M{z)  = I.  Thus  if  M is  invertible,  then  for  any  stable 
basic  set  (p  for  V,  there  is  a unique  dual  basic  set  ip  for  W.  □ 

We  remark  that  from  the  definition,  any  biorthogonal  basic  sets  must 
have  linearly  independent  integer  translates.  Thus  if  M as  in  Theorem  2.1  is 
invertible,  any  basic  set  for  V or  W which  is  stable  must  in  fact  have  linearly 
independent  integer  translates. 

Now  for  m E 7Z,  m > 2,  we  say  a space  V of  functions  on  1R  is  m-refinable 
if 

f £V  =>/{-)£  V. 

m 


M(z) 


E aM*)** 


Lfc=— oo 
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If  <j>  is  a basic  set  for  an  m-refinable  space  V,  then  is  in  V,  and  so 

OO 

^(-)  = E Ok<K.-k) 

TTI 

K — — OO 

for  some  finite  set  of  matrices  C*,.  Thus,  <j>  satisfies  the  refinement  equation 

OO 

^ = E - k). 

k=—oo 

In  this  paper  we  shall,  for  simplicity,  consider  only  the  case  m = 2. 

§3.  Biorthogonal  Refinable  Splines 

For  any  integer  n > 1 we  let  Nn  denote  the  uniform  B-spline  of  degree  n with 
simple  knots  at  0, 1, . . . , n + 1.  We  now  fix  n,  and  define 

= Nn( 3x),  </>2(a;)  = Nn(3x  - 1),  <f>3(x)  = Nn(3x  - 2). 

Then  {4>i , fa , <j>3  } is  a stable  basic  set  for  the  space  V of  all  spline  functions 
of  degree  n with  compact  support  and  simple  knots  at  \TL-  Clearly  V is 
refinable.  We  wish  to  find  a refinable  space  W of  spline  functions  of  compact 
support  which  has  a basic  set  which  is  dual  to  {<f>i,<j>2,4>3}-  From  Theorem  2 
we  see  that  this  is  equivalent  to  finding  a basic  set  {V’i>  1P3}  for  W such 

that  the  matrix  M in  (6)  is  invertible. 

We  shall  choose 

■^i(s)  = Nn( 2x),  ip2(x)  = Nn(2x  - 1), 

and  V>3  to  be  a spline  function  of  degree  n with  knots  in  and  support  in 
[| , 2n  — |] . Thus  W is  a space  of  spline  functions  of  degree  n and  simple  knots 
in  j7Z,  which  contains  all  spline  functions  of  degree  n with  compact  support 
and  simple  knots  in  |2Z.  For  any  function  / in  W,  /(|)  has  knots  in  \ TL  and 
so  lies  in  W.  Thus  W is  refinable. 

Theorem  3.  We  can  choose  as  above  so  that  there  are  biorthogonal  basic 
sets  for  V and  W , or  equivalently  that  there  is  a unique  basic  set  for  W dual 
to  the  basic  set 

Proof:  The  space  Wo  of  spline  functions  of  degree  n with  simple  knots  in  |2Z 
and  support  in  [0, 2 n]  has  dimension  3n,  while  the  space  W\  of  spline  functions 
of  degree  n with  simple  knots  in  |ZZ  and  support  in  [0, 2n]  has  dimension  In. 
Thus  we  may  choose  linearly  independent  functions  fi,  ■ ■ ■ , fin  in  Wl  with 
support  in  [i,  2 n — |]  which  together  with  Wq  span  Wi . We  write 


4 n 

i>3  = E°*/*- 

fc=i 
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It  remains  to  choose  ai,...,ain  so  that  M in  (6)  is  invertible.  Now 
Mij(k)  = 0 except  in  the  following  cases:  Mn(k),  — ^ < k < M]2(T), 

_n±i  < k < a=i;  M13(k),  -2n  + l < k < §;  M21(k),  -f  < k < 
M22(fc),  < fc  < f;  M23(fc),  -2n+l  < fc  < *±1;  M31(fc),  -»=*  < fc  < 

“f2;  M32(fc),  -f  < k < JSfi;  M33(k),  -2n  + 1 < fc  < ^.  Then 

n 

det  M(z)  = X bkZk, 
fc=- 3n+l 

for  some  numbers  6^ , — 3n  + 1 < k < n,  which  are  linear  functions  of 
oi,...,a4„.  The  condition  det  M(z)  = 1 then  gives  4n  linear  equations  in 
4n  unknowns  <xj , . . . , a4n  and  we  shall  show  that  this  system  is  non-singular. 

Suppose,  to  the  contrary,  that  the  system  is  singular.  Then  we  may 
choose  ai, . . . , a4„,  not  all  zero,  so  that  det  M(z)  = 0.  Then  the  columns  of 
M are  linearly  dependent  in  the  sense  that  there  are  Laurent  polynomials  pi , 
P2i  P3,  not  all  zero,  so  that 

3 

X M(z)ijPj(z)  = 0,  i = 1, 2, 3. 

j = 1 

Writing 

oo 

Pj{z)  = X ci(k)zk’  J = 1-2>3> 

k=  — oo 

this  becomes 

3 oo  oo 

X X M^k)zk  X c^)z‘=  0,  ’ = 1-2-3’ 


XX  Mij(l)cj(k  — l)  = 0,  * = 1,2, 3,  k € ZZ, 

J=1  i=— oo 

which  on  recalling  (5)  gives 


/: 


°,  * = 1,2,3, 


k e7Z, 


/ = X X 

j=l  l=—oc 

Thus  / is  orthogonal  to  V.  Note  that  since  the  integer  translates  of 
'01,  02,  03  are  linearly  independent,  / is  not  identically  zero.  Since  only  a 
finite  number  of  coefficients  j = 1,2,3,  l € ZZ,  are  non-zero,  / has 
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compact  support.  Suppose  that  the  support  is  [a,0  + 2 n],  a,  0 £ 7Z,  but  not 
in  [a  + 1,0  + 2n]  or  [a,  0 + 2n  — 1],  It  is  easily  seen  that  a < 0 and 

0 

f = ~ j)  + 9, 

j=a 

where  ca  ^ 0 ± cp  and  g has  support  in  [a,  0 + 2n]  and  lies  in  the  space  of 
splines  of  degree  n with  knots  in  |2Z,  which  we  shall  denote  by  Z. 

We  first  note  that  for  some  hi  in  Z with  support  in  [0  + ^ — 1,0  + 2 n], 

m :=  cpijj3(.  -0)  + hx 

is  orthogonal  to  those  elements  of  V with  support  in  [0  + 2n  — 1,  oo).  Now 

1 /3_1 
& j—a 

is  orthogonal  to  those  elements  of  V with  support  in  [0  + 2 n — 2,  oo)  and  on 
this  interval  coincides  with  a function 


V2  ■=  c0ip3(.  - 0)  + h2, 

where  h2  is  in  Z with  support  in  \0  + ~ — 2, 0 + 2n],  Continuing  in  this  way 
we  recursively  construct 


Vj  ■=  cpip3(.  - 0)  + hj,  j = 1, . . . , 4n, 

which  is  orthogonal  to  those  elements  of  V with  support  in  [0  + 2 n — j,  oo), 
where  hj  is  in  Z with  support  in  [0  + ^ - j,0  + 2 n}.  In  particular,  r)in  has 
support  in  [0  — ^r,0  + 2n]  and  is  orthogonal  to  those  elements  of  V with 
support  in  [0  — 2n,  oo). 

Choose  F with  _F(n+1)  = ?j4n  and  with  support  in  (— oo, 0 + 2n].  Now  for 
j = 0, . . . , 12n  — 1,  let  Bj  be  the  B-spline  of  degree  n with  knots  0 — 2n  + 

0 — 2 n+  ...,/?  — 2 n+  ^+g+1.  Then  since  Bj  is  in  V, 

/OO  rOO 

BjTHn  = / BjFn+1 

-OO  J — OO 

= [0-2n+3-,0-2n  + J-±±,...,0-2n+3+n3  + 1]F. 

Thus  F vanishes  at  0 — 2n  + j — 0, . . . , 12n  — 1.  Now  F coincides  on 
[0  — 2n,0+2n]  with  a spline  G of  degree  2n+l  with  support  [0—3n—^,0+2n] 
with  knots  at 

0 — 3 n ~ 2’@  — ® n ’ ■ ■ ■ i0  ~ 2’^'^  + 4’^  + 2'"’'^  Jr^‘n  ~ 4 ’0  + <3,n' 
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It  then  follows  from  the  Schoenberg- Whitney  Theorem  [5]  that  G vanishes 
identically  on  [/?  — 2n,  (3  + 2 n]  and  hence  so  does  7/4 n.  So  774  n has  support  in 
\P  — — 2 n]  with  knots  in  \TL,  and  so  774,,  vanishes  identically,  which  is 

a contradiction.  Thus  the  linear  system  is  non-singular,  which  completes  the 
proof.  □ 

Finally  we  note  that  if  1P3  is  as  in  Theorem  3 and  M is  given  by  (6),  then 
the  basic  set  for  W dual  to  {4>\ , 4>:i}  is  {Vh, ^2, ^3}  given  by 

OO 

fc=— 00 


where 


OO 

Y,  C^k  = (W1)*  = (adJ  = (adJ  M(z-1))T. 

k=— 00 
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